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We find the first example of a quantum Berenzinskii-Kosterlitz-Thouless (BKT) phase transition 
in two spatial dimensions via holography. This transition occurs in the D3/D5 system at nonzero 
density and magnetic field. At any nonzero temperature, the BKT scaling is destroyed and the 
transition becomes second order with mean-field exponents. We go on to conjecture about the 
generality of quantum BKT transitions in two spatial dimensions. 
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Introduction. — Holography [IH3] has become an impor- 
tant tool in the investigation of strongly-coupled systems. 
Despite being restricted to a special class of theories — 
those with simple gravitational duals — the technique has 
found interesting applications to the physics of the quark- 
gluon plasma [3]. Recently, there have been attempts to 
use holographic models to approach condensed matter 
systems — like the Fermi gas at unitarity [5, 6j, superflu- 
ids [7] and non-Fermi liquids [HHTT] — with various degrees 
of success. 

Phase transitions are one central physical concept that 
can be studied holographically. They are present in many 
holographic models and frequently allow simple geomet- 
ric interpretations. However, most holographic phase 
transitions are either hrst-order or second-order 

with mean-field exponents. The reason for the mean- 
field behavior is the large N parameter which suppresses 
quantum fluctuations in the gravity theory and allows the 
latter to be treated semiclassically. On the other hand, 
the most interesting questions are usually beyond the 
mean- field approximation |15| . or beyond the Landau- 
Ginzburg- Wilson paradigm altogether |16j . It is very im- 
portant to investigate whether the holographic method 
can be used to study these non-mean-field phase transi- 
tions. 

In this Letter, we show that, even within the confines 
of large-iV field theories, another type of phase transi- 
tion is possible — namely, those with the scaling behavior 
of the Bcrczinskii-Kostcrlitz-Thouless (BKT) phase tran- 
sition pTl [18] . Recall that in the BKT phase transition, 
the order parameter scales as exp(— cj\JT c — T) near the 
critical temperature T c 19 . In our case, interestingly, 
the phase transition is a quantum phase transition, occur- 
ring at zero temperature in 2+1 spacetime dimensions. 
The explicit example we consider is a 2+1 dimensional 
theory at a finite density d of a conserved charge and 
a magnetic field B. At a particular value of the "fill- 
ing fraction" v = d/B, the system suffers a transition 
to a broken symmetry state, and the scaling of the or- 
der parameter is the same as in a BKT phase transition, 
a ~ exp(— c/yfv c — v). We shall call this phase transi- 
tion the "holographic BKT phase transition," although 



it happens in a context different from the original BKT 
phase transition. The BKT scaling occurs in quantum 
mechanics with a 1 /r 2 potential, and has been speculated 
to describe the chiral phase transition in QCD with large 
number of flavors [20) . 

On the gravitational side, the transition occurs due to 
the violation of the Breitcnlohner-Freedman bound [2T] 
in the infrared region by the scalar field dual to the or- 
der parameter |22). Although this mechanism seems sim- 
ilar to that of a second-order Landau phase transition, 
the field in the bulk actually represents an infinite tower 
of states in the boundary quantum field theory. At the 
phase transition, an infinite number of field theory modes 
become unstable at the same time — an extremely unnat- 
ural situation within the Landau theory. 

The gravitational description of the holographic BKT 
transition is given in terms of a probe brane minimizing 
its worldvolume action in a fixed geometry. In Ref. [23], 
three of us studied one such probe describing a 3+1 di- 
mensional field theory. In 3+1 dimensions the competi- 
tion between finite density and the magnetic field gave 
rise to a second-order transition with mean-field scaling. 
(See Ref. [33] HI] f° r technical details and a discussion of 
the related literature.) In this work, we will show that the 
same setup in 2+1 dimensions gives rise to BKT scaling. 
The fact that the density and the magnetic field have the 
same mass dimension in 2+1 dimensions will be crucial. 

The D3/D5 system.— We write the AdS 5 x S 5 back- 
ground metric as 

9=(r 2 + y 2 )(-dt 2 + dx 2 ) 

+ ^—(dr 2 +r 2 dn 2 + dy 2 +y 2 dn 2 ), (1) 
r z + y z 

where a? is a spatial three-vector and dQ?. is the met- 
ric of a unit-radius two-sphere. In these coordinates the 
boundary of the AdSs is located at r 2 + y 2 — » oo and the 
horizon at y = r = 0. 

We consider Nf probe D5 branes in this geometry, 
making an Ansatz that the probes wrap the first two- 
sphere, the "radial coordinate" r, time, and two spatial 
directions. The branes have a profile parameterized by 



2 



V = y{ r )- The dual field theory is Af = 4 SU(N) super- 
Yang-Mills theory at large 't Hooft coupling A coupled 
to Nf fundamental hypermultiplets along a 2+1 dimen- 
sional defect. This theory has a 17(1) x SU(2) 1 x SU{2) 2 
global symmetry, where the two SU(2) factors are chi- 
ral R-symmetries. The field y is dual to a condensate of 
the field theory which can spontaneously break the sec- 
ond SU(2) factor. We are interested in the chiral phase 
transition between y = and y ^ configurations. 

We turn on a density of the U(l) flavor charge and 
a magnetic field coupled to this charge. On the gravity 
side, this corresponds to turning on a field strength for 
the U(l) gauge field dual to the current, F = A' t {r)dt A 
dr + Bdx 1 Adx 2 . The radial electric field A' t is supported 
by charge behind the AdS horizon, so these branes must 
extend to the bottom of AdS. 

The Dirac-Born-Infcld (DBI) action density for these 
probe branes takes the form 



behavior near the phase transition. In this regime, the 
embedding can be found by matching the solutions in 
two regions: a small-r nonlinear core and a large-r linear 
tail. 

In the core region, we can set B 2 = p 2 /7, and neglect 
r 4 compared to p 2 . The action then becomes 



S 



-Afp I dr VTT y 




7(r 2 + y 2 ) 2 



(6) 



Our brane embeddings, which extremize this action and 
obey the boundary condition y(0) = 0, form a one- 
parameter family of solutions related by scaling: y^(r) = 
£/( r /£); where y = f(r) is one particular embedding. 
At large r, any solution in the family goes to the linear 
regime, where 



f(r) = Vf(-a + ailog r), r > 1. 



(7) 



S = -N j drr 2 ^Jl + y' 2 - A' 2 J I 



B 2 



(r 2 +y 2 ) 



2\2 



(2) 



where Af = ^N f N c /(2ir 3 ). The field A t only appears 
through derivatives, so d = SS/SA' t is r-indepcndcnt. It 
is the density of the U(l) charge. The action of y at 
fixed density is obtained by Legendre transforming |2]) 
with respect to A' t , 



i / Af 2 r i B 2 

S = I dr d 2 + Af 2 r 4 + — . (3) 



(r 2 + y 2 ) 



The onset of the phase transition can be found by ana- 
lyzing the stability of small perturbations around y = 0, 
which are described by the quadratic part of 



2 r 2 y / p 2 +B 2 



(4) 



where p = d/Af. We pause to note two features of this 
Lagrangian. Near the boundary at large r, y/r fluctu- 
ates as a scalar with mass squared m 2 — — 2 in AdS.4. 
However, at small r (the IR of the dual theory), y/r be- 
haves like a scalar in AdS 2 with m 2 = —2B 2 /(B 2 + p 2 ). 
If p/B < \/7, then the mass of y/r in the IR region is 
below the Breitenlohner-Freedman bound of stability for 
the effective AdS2, ™| F = — 1/4. In this case the trivial 
embedding y = is unstable and the ground state should 
instead have y 7^ 0. Thus, the chiral phase transition oc- 
curs at the filling fraction 



_d_ 



7\N f N c 
~2^3 



(5) 



Numerically solving the equation of motion by shooting, 
we find a solution / with ao = —.211 and a\ — .0585. A 
general embedding in the core then has the asymptotics 



To 



r » r = £e 



(8) 



It would be incorrect to use this expression for the linear 
tail too far away from the core: the terms neglected when 
the action ^ is written down become important. But 
in the linear regime, the two independent solutions for y 
can be found without additional approximation, 



f ± {u)=u^' 2 2 F 1 
where 



1 ± ia 3 ± ia ia 
~ 8" '~ s - ' T'" 



VpI - p 2 

Vp 2 + B 2 



( p 2 +B 2)l/4 



(9) 



(10) 



When the filling fraction is just below the phase transi- 
tion, a is real and small. From the two solutions found in 
Eq. (9), one can construct a particular linear combination 
which asymptotes to 1/u near the boundary, 



U(u) = c + f + (u) + c_/_(u) 



1 

u 



(11) 



This choice amounts to choosing zero bare mass for the 
dual flavor. The coefficient of the 1/u term is propor- 
tional to the condensate of the dual theory. Denoting 
the condensate as a, the solution in the linear regime is 
y(r) = -(p 2 + B a )-VVn(«). 

For small u and small a, /„ can be expanded 



In the broken phase, the condensate can be found by 
solving the equation for the embedding with appropriate 
boundary conditions. We are interested in the critical 
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(l+ia)/2 



(12) 
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FIG. 1. A plot of the condensate as a function of magnetic 
field at zero and finite temperature near the zero-temperature 
transition. The dashed black line indicates zero-temperature 



numerical data and the solid blue line our prediction, Eq. (171, 
computed further to C(a 3 ). The color dashed curves repre- 
sent numerical data at temperatures of T = — p 1 ^ 2 x lO^ 11 
(left) and T = |p 1/2 x 10" 10 (right). At any nonzero tem- 
perature, the condensate scales with a mean-field exponent 
near the transition and then asymptotes to the BKT scaling 
at large magnetic field. 



where 



Ci = 



r 2 (i/4) 

23/4^3/2 



2V 4 r 2 (5/4)(ln256 - 7 r) 

62 - ZS/2 ■ i 13 ) 



To leading order in a, we rewrite Eq. (12 1 as 



y(r) 



2C lP 



(P 2 



n = (p 



B2) 3/8-- Sln U lE - 



2 ^l/4 I 2C 2 

B J exp 



(14) 
(15) 



To match the core and linear solutions at r = £, the 



argument of the sin in Eq. ( 14 ) has to make 7r between 
r o ~ £ an d fi ~ s/B. To exponential accuracy, we then 
have £ ~ e~ 2 ^^ a . Comparing the prefactor of the expo- 
nent, we find 



(16) 



Thus, the phase transition in the D3/D5 system obeys 
BKT scaling. The exponent as well as the prefactors can 
also be found by matching both y and y' at £, 



B 2 exp 



9i 



2ai 



0(a) 



(17) 

We also compare this result to numerical data in Fig. [I] 
The way we matched the core and the tail of the em- 
bedding makes it clear that there exists an infinite set of 
"Efimov extrema," which we name for their resemblance 
to Efimov states [25]. Indeed, the most general match has 



an argument nir in the sin of Eq. ( 14 ) between r ~ ro 
and r ~ n, for n a positive integer. These embeddings 
go through n/2 oscillations between the bottom of the 
brane and the boundary. Matching y and y' at r = £ 
as before, we find a condensate corresponding to these 
extrema of a n ~ e - ""/", where <7i is the condensate for 
the simplest embedding above. The additional extrema 
are not ground states, however: if one denotes by Fq the 
free energy of the trivial embedding, the free energy of 
the n-th extremum goes like F — F n ~ e -2nn/a^ g0 ^. ue 
n = 1 embedding is the ground state. The infinite tower 
of extrema, spaced by the same factor exp(— 2n/a), is 
reminiscent of the Efimov states. 

We next examine the D3/D5 system at finite tem- 
perature T. To do so, we embed the D5 probes in an 
AdS-Schwarzschild black hole geometry. The embedding 
equations can now only be solved numerically. However, 
at any temperature the magnetic field does not contribute 
in the effective AdS2 region in the IR. We could there- 
fore predict that the BKT scaling is lost for T > 0. We 
numerically solved the embedding equations at a num- 
ber of different temperatures. (The techniques employed 
are similar to those used in Ref. [23] .) We plot some of 
our results for the condensate in Fig. [l] At fixed finite 
temperature, we find that there is still a chiral symme- 
try breaking transition but that it is a mean-field second 
order one. Far away from the transition, we recover the 
exponential BKT scaling at zero temperature. 

Lifshitz systems. — In the example that has been just 
considered, the BKT scaling arises from an IR AdS2 re- 
gion and the existence of a scalar with mass squared 
crossing the BF bound at the transition. We may ex- 
pect the BKT scaling to be present in many other cases. 
We now demonstrate that the BKT scaling indeed con- 
trols the critical behavior in phase transitions in a rather 
large class of models. 

Our models are built on the so-called Lifshitz geome- 
tries |26j . which are invariant under anisotropic scaling, 
t — > X z t } x — > Ax. Under holography, a geometry of this 
type is conjectured to be dual to a scale-invariant bound- 
ary field theory with dynamic critical exponent z. Con- 
structing explicit Lifshitz solutions from string theory is 
quite involved llj; our approach here is phenomenologi- 
cal. Nevertheless, once a few conditions (to be specified 
below) are satisfied, the BKT scaling is rather generic — 
which seems to indicate that it should occur in string- 
theoretical realizations of Lifshitz geometries as well. 

We consider probe branes with an induced metric 



p [ff] = -(r 2 +y 2 ) z dt 2 +(r 2 +y 2 )dx 2 + 



dr 2 (l + y' 2 ) + r 2 g int 
t 2 + y 2 

and constant dilaton. Here, r is the "radial coordinate" 
of the geometry, z is the dynamical exponent of the dual 
theory, x G BP is a spatial vector, <7i nt is the metric on a 
wrapped space of dimension fc, and y = y(r) is the profile 
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of the branes. The D3/D5 system corresponds to z = 1, 
D = 2, and k = 2. 

We now turn on a density and magnetic field in the 
field theory, dual to U(l) field strengths on the branes. 
Assuming that Chern-Simons terms play no role, the ac- 
tion of the probes at fixed density is 



S=-Af z drR z ^/l + y' 2 




r 2k R 2(D-k) 1 



(19) 



r ' 2 + y 2 , M z is a prefactor, 



where we have defined R 
and p is a rescaled density d = M z p. 

As before, we attempt to find the location of the phase 
transition by looking at small fluctuations around the 
trivial embedding y(r) = 0. In the UV, y/r is a massive 
scalar in AdSz+3, corresponding to an operator with di- 
mension A = z/2 + y/(z + 2) 2 - 4fe/2 + 1 in field theory. 
At small r, y fluctuates in two different ways depending 
on the value of D. For two spatial dimensions, y/r fluctu- 
ates as a m 2 = z—l — kB 2 /(B 2 +p 2 ) scalar in AdS z+ i; for 
D > 2, y/r behaves like a m 2 = z — 1 scalar in AdS z +i. 
The holographic BKT transition is then only possible for 
D = 2, where a density and magnetic field have the same 
dimension. The transition occurs at a critical density of 



P 



(20) 



for which the bound 4fc > z 2 must also be satisfied. This 
is also the bound for which a magnetic field at zero den- 
sity will break chiral symmetry. 

Using the same method as employed in the D3/D5 
case, we study zero bare mass by taking y to be normal- 
izable at large r with y ~ 1/r . Holographic renor- 
malization relates the coefficient of this term to the con- 
densate er, which we find to be 



( i o 2 + B 2 )T(A-l)exp 



u zCi zao \ 
(21) 

where Ci, C2 are complicated constants that depend on 
A and z, and oq, a\ are the asymptotic data of a core 
solution to which we match. We can therefore find BKT 
scaling for any dynamical exponent, provided that the 
magnetic field can break chiral symmetry. 

Discussion. — We therefore see that holographic phase 
transitions in 2+1 dimensions, with the ratio of the mag- 
netic field and the density as a control parameter, may 
generally be of the BKT type. It would be extremely 
interesting to find such a BKT transition away from the 
large N and strong coupling limits. 

The existence of an infinite number of "Efimov vacua" 
is clearly related to the breakdown of the Landau's the- 
ory. In Landau's theory of phase transitions, the order 
parameter is the only low-energy degree of freedom. In 



our case, the role of the order parameter is played by 
the whole function y(r), where r is usually interpreted 
as the energy scale. One needs to futher investigate the 
low-energy effective theory of the D3/D5 system. 
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